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Abstract 

We analyze the relation between the concept of auxiliary variables 
and the Inverse problem of the calculus of variations to construct a 
Lagrangian from a given set of equations of motion. The problem of 
the construction of a consistent second order dynamics from a given 
first order dynamics is investigated. At the level of equations of motion 
we find that this reduction process is consistent provided that the 
mapping of the boundary data be taken properly into account. At 
the level of the variational principle we analyze the obstructions to 
construct a second order Lagrangian from a first order one and give 
an explicit formal non-local Lagrangian that reproduce the second 
order projected dynamics. Finally we apply our ideas to the so called 
"Noncommutative classical dynamics". 



1 Introduction 

Inspired from recent results in string theory some authors, have 
proposed a "noncommutative quantum mechanics on the plane" j2j. 
Here by quantum mechanics they mean a classical field theory for 
the classical "probability field" ip{x) defined in a space where the 
coordinates do not commute 

[x i ,x j ]=i6 ij , (1) 
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where 9 lJ is a constant antisymmetric matrix. The tools used recently 
for the analysis of noncommutative field theory can be applied to this 
case. In particular, the introduction of noncommutativity in field the- 
ories renders the classical field functions on space-time to operators or 
matrices whose treatment resemble very nearly the techniques usually 
employed to quantize classical field theories [3]. It is also possible to 
use the Weyl or Moyal star product in a constructive way to "deform" 
the standard product of functions. Indeed, a procedure to obtain, 
in the appropriate limit, this non-commutative quantum mechanics 
starting from a noncommutative field theory was given by 4 . The 
complete picture that emerges from here is quite non-local and has 
been investigated only for some simple cases jHj. 

Another beautiful approach to obtain a consistent quantum dy- 
namics in phase space with a curved symplectic structure, in contrast 
to the noncommutative space-time (0), is deformation quantization 
0. Here we start with an associative commutative algebra A and de- 
form it into a new associative but noncommutative new algebra. This 
can be done using H as a deformation parameter but also we can use 
a 13 as a deformation parameter (see for example [B]) where 

{:'.: J } (2) 

Here {z l , z 3 } denotes the Poisson bracket in the "phase space" defined 
by the coordinates z l and a %3 is a constant matrix. In this case, we can 
use an ordering prescription, say symmetric Weyl ordering, to get a 
consistent picture of the algebra of observables by replacing the stan- 
dard product of operators by the star or Moyal product in the space 
of its associated symbols 9 a . The picture that emerges from here is a 
quantum mechanics on phase space defined by some Hamiltonian and 
the structure given by (|2j . Of course this approach is very different 
from the one mentioned above. To recover from deformation quanti- 
zation a Schroedinger equation and a wave function in configuration 
space, when the coordinates do not commute, is not an easy task. The 
reason is that the usual momenta associated with the noncommutative 
coordinates are not auxiliary variables, in the usual sense. The conse- 
quences that this fact has upon the resulting quantum mechanics are 
quite severe. Indeed, it is not known how to project, in a consistent 
way, the Wigner distribution functions associated with a given prob- 
lem to distributions of probability in configuration space. For non 
constant a 13 the analysis is much more involved. Very interensting 
recent results along this line started with the seminal work jlOj and 
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have been applied to a sigma model in 

One obvious question is if there exists a noncommutative classi- 
cal mechanics that upon some quantization procedure, say canonical 
quantization, gives the "noncommutative quantum mechanics on the 
plane". Since Classical Mechanics (in its first order formulation) is 
covariant with respect to a wide class of transformations -Darboux 
transformations- the question of noncommutativity in Classical Me- 
chanics does not make sense unless we fix coordinates using some 
physical criteria. 

An interesting approach not pursued here could be to consider 
classical dynamics with a noncommutative time. The result of such 
analysis could be similar to the non-local classical dynamics analyzed 
in P2|. 

Here we will consider a different problem which is not related, at 
least in any direct way, to the approach of constructing noncommu- 
tative field theories 1 . The approach that we will present here could 
be related with deformation quantization. The first problem, at a 
Classical level, is how to define in a consistent way the configuration 
space dynamics when the coordinates on this space do not commute 
in the symplectic sense (J2J). Namely, to define from a given first order 
variational principle the associated second order variational principle 
in configuration space upon the elimination of momenta from the first 
order dynamics. As far as we know, there is not a complete answer 
to this problem. A previous analysis of the Inverse problem from this 
perspective can be found in [22] . 

In its simplest setting the problem can be stated as follows: giving 
a classical system described by some Hamiltonian of the form H = 
T + V, and some (constant) symplectic structure, can we define the 
second order equations of motion in configuration space in such a way 
that to every solution of the original first order system of differential 
equations we can associate a solution of the second order equations 
in configuration space? And, as a second step, can we define the 
associated quantum theory? 

A previous attempt to analyze the problem of the quantization of 
classical systems with nonstandard symplectic structures can be found 
in JH]) where the importance of the existence of a second order La- 
grangian for the quantization is remarked. In this context a discusion 
about the restricted inverse problem of the calculus of variations and 

1 A recent attempt to construct a noncommutative field theory with fields that are 
noncommuting can be found in [H] 
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the quatntization of first order systems can be found in ^7] . 

The aim of this note is to give an answer to the first of these 
questions and to show that there are still some uncovered problems. 
The analysis of the second question will be presented elsewhere. 

For that end we will review in section 2 the basic ideas of aux- 
iliary variables in first order systems and the inverse problem of the 
calculus of variations, i.e, the necessary and sufficient conditions for 
the existence of a Lagrangian for a given set of first order dynamical 
equations of motion ^3]. This problem is of relevance because there 
is a relation between the existence of auxiliary variables in the first 
order dynamics and the existence of a Lagrangian function for the 
configuration space equations of motion. Here we want to confront 
the following paradox: From one hand a Lagrangian for an integrable 
system of first order differential equations always exist |13| . From the 
other hand we can always associate to a given system of second order 
differential equations a system of first order differential equations. By 
eliminating the auxiliary variables from the first order system we can 
recover the second order system and by the same procedure obtain 
also a Lagrangian for these second order equations, using the first or- 
der variational principle. But it is well known that the Lagrangian for 
the second order differential equations does not always exist ^j. So, 
which are the general conditions over the first order Lagrangian func- 
tion to give, upon reduction, a second order Lagrangian associated to 
an equivalent second order system? As we will see, this problem can 
be formally solved with a non-local Lagrangian that gives an equiva- 
lent set of equations of motion in the reduced space. The Lagrangian 
equations of motion in the reduced space are equivalent to a set of 
second order equations of motion, but the equivalence with the orig- 
inal first order dynamics can be implemented only through a proper 
consideration of the boundary data and its "projection" to the second 
order dynamics. 

Given the coordinates of the configuration space we will investigate 
in section 3 the case when we have at hand a Hamiltonian to describe 
the first order dynamics. In section 4 we will present a procedure to 
construct a non-local Lagrangian for theories with a given Hamiltonian 
of the form H = T + V and a constant symplectic structure. In 
section 5 we will define the "noncommutative classical dynamics" . We 
will analyze the precise equivalence between "Newtonian equations of 
motion" and the first order formulation of the theory. Section 6 is 
devoted to examples and section 7 to conlusions and perspectives. 
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2 First order Lagrangians and the in- 
verse problem of the calculus of varia- 
tions 

In this section we will review two central ideas relevant for our discus- 
sion. The first one is the definition of auxiliary variables for a given 
variational principle .18; • The second one is about the existence of a 
Lagrangian for a given first order system of differential equations ^3] • 
Then we will relate this two ideas to construct a Lagrangian for the 
reduced system that results from the elimination of the auxiliary vari- 
ables in the original system and point out the possible obstructions to 
perform such procedure. 



2.1 Auxiliary variables 

Consider the action 

S= dtL(x, x, x, ....), (3) 



where the Lagrangian depends on the coordinates x a ,a = 1,2...N 
and its derivatives with respect to time up to any finite order. The 
equations of motions will be denoted by 

n=0 

where D is the total derivative with respect to time given by 
^ . d .. d d 

D = x Tx +x Tx + - + Wv (5) 

and x( n ) denotes the re— derivative with respect to time of the coor- 
dinate x. Now suppose that the coordinates can be splited into two 
groups 

x a ^(y\z a ), (6) 

such that 

^=0^z a =Z a (y\y\y\...). (7) 

Here we want to stress that the procedure to obtain the variables z a 
in terms of the reduced space variables y % is purely algebraic 2 . In 



2 There are some particular examples where this step can be realized by integration. As 
far as we know there is no a systematic generalization of the ideas below for this case. 
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that case the variables z a are called auxiliary variables. The reduced 
Lagrangian L R is defined by 

L R (y\y\y\....)=L(y\Z a ,y\Z a , ....). (8) 

Then we can state the following Theorem: The original equations of 
motion 

5L - 5L - (9) 

5y l 5z a 

are equivalent to the equations 

^ = 0, z a = Z a . (10) 

The proof is straightforward using the chain rule 
5L R _ 5L 5L dZ a d / 6L dZ a \ k d( fc ) ( 5L dZ a \ 



<fy* Sy i 5z a dy i dt \5z a dy 1 J " di( fc ) \Sz a Q y {ky J ' 

(11) 

The equations = imply z a = Z a (by asumption) and = 

(by pij))- Conversely, the equations = and z Q = Z a imply 

^ = (by assumption) and 44 = fbv (jTT|> ). 

As an example lets see how this works with the standard Hamil- 
tonian dynamics. The first order Lagrangian is 

L = xp — -p 2 — V(x), 

while the reduced Lagrangian (after the elimination of the momenta) 
is 

L R = h 2 -V{x). 
To see the equivalence we write the chain rule 

x + V'{x) = — -- — , 
ox dt dp 

which shows explicitly the required equivalence. It is interesting to 
notice that if we fix as initial conditions x,p at some initial time to in 
the equations of motion associated to L(x,p) then the initial conditions 
for the reduced system 'x + V'{x) = are the standard Newtonian 
position and velocity at initial time to. 
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2.2 Inverse problem for the calculus of varia- 
tions for first order systems and reduction 

In this subsection we will construct a criteria for the existence of a 
second order Lagrangian given a first order one based on the concept 
of auxiliary variables. 

By splitting the original variables x a into z a and y % and by choosing 
the y l as the variables that define the reduced space we can write a 
relation of the form 

where Gi are linear combinations of the first order equations of motion 
associated to the original first order Lagrangian L and its derivatives. 
The question is if this linear combinations can be written as a vari- 
ational derivatives with respect to the reduced variables y l for some 
function Lft(y, y), i.e., if the equations of motion of the reduced space 
Mi(y,y,y) = comes from a variational principle. This will be true 
iff 

SMjjt) _ 5M 3 (t>) 

8yi(t>) 5f(t) ' 1 ' 

These conditions are the Helmholtz conditions for the restricted in- 
verse problem of the calculus of variations associated to the system 
Mi = 0. If this conditions are meet then Gi are variational total 
derivatives. A generalization of this problem is the question about 
the existence of a matrix A i3 (y, y) with det A ^ such that the equiv- 
alent system A^Mj = can be deduced from a variational principle, 
i.e., if an "integrating factor" exists such that the relations (1121) are 
total variational derivatives. This is the unrestricted inverse problem 
of the calculus of variation. A review about this interesting problem 
can be found in |14j . 

The conditions l)13j) can be written in terms of the original first 
order Lagrangian L through the functions Gi to see which are the 
conditions under the Lagrangian L to give a second order variational 
formulation of the system Mi upon the reduction process. 

If the variables used in the reduction process are auxiliary then Gi 
is a total functional derivate and the conditions (jl<-i|) for the existence 
of a Lagrangian are automatically satisfied. The reduced Lagrangian 
Lr is a Lagrangian for the reduced second order equations of motion. 

The functions Gi contains the precise information about the alge- 
braic manipulations performed on the original system of equations to 
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arrive at the second order ones in the reduced space. The relations 
(fT2|) do not imply that the two sets of equations |^ = 0, j^t = 

and Mi = 0, Mk = are equivalent. The equivalence between these 
two sets of equations depends crucially on the structure of the orig- 
inal equations of motion and on the reduction process. There are 
two interesting scenarios: a) reduction with auxiliary variables while 
maintaining the equivalence (as in the standard Hamilton to Lagrange 
description in Classical Mechanics), b) reduction with non-auxiliary 
variables while maintaining the equivalence (a case where the second 
order Lagrangian does not exist (in the restricted sense stated above)). 

If the variables used in the reduction are not auxiliary Gi is not 
a total variational derivative and the Lagrangian associated to the 
system M, does not exist. But it is still possible that a Lagrangian as- 
sociated with the equivalent system A' J Mj = could be constructed. 
We will not pursue this general problem here but we will take a short- 
cut to construct a non-local Lagrangian in the case where we do not 
have auxiliary variables to perform the reduction procedure and the 
symplectic structure is constant. 

The analysis of the unrestricted inverse problem of the calculus 
of variations for first order system is, in contrast, much more easy to 
workout. It can be stated as follows |13j : Given a first order system 
of differential equations in a space defined by the coordinates x a ,a = 
l,2...2n 

x a - f a (x) = 0, (14) 

we can construct a Lagrangian for this given set of equations in the 
form 

^- = a ab (x)(± b -f b ), (15) 

where a a i, is a non singular antisymmetric matrix and deter ^ 0. This 
Lagrangian can be constructed by the following procedure [Ej. First 
write a general Lagrangian function as 

L = t a (x)(x a -r). (16) 

Now, using (fT3|) we see that £ a (x) must satisfy 

where Cf is the Lie derivative along the vector field f a and the matrix 
a is given by 

__ di b d£ a 
aab ~^~ (18) 
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in order that the equations of motion associated with the Lagrangian 
((TB|) be equivalent to the set of the given equations Q14j>. The condition 
l|17jl always has solution (at least locally) 13; so it is always possible to 
construct a variational principle for a given set of first order differential 
equations. 

In particular, a Hamiltonian form for the given system (|14|) can be 
constructed if we can find a function H{x) such that 

B ff 

£ a f a = H(x), a ab f b = —. (19) 

Conversely, given H(x) satisfaying (fUJ)) the function £ a satisfy (|17|). 
The symplectic structure of the resulting Hamiltonian system is given 
by a ab {a = dt). 

We can always choose coordinates (at least locally) such that the 
given system (|14[) can be transformed to the "normal form" 

± a - f a ( x ) = o c - e +n = o, e +n - f +n (o = o. (20) 

Using these variables and taking as the reduced space we obtain, 
upon reduction, the second order system 

e-r(e,i 3 ) = o, F i (e,e) = f+ n (o ... (21) 

The Lagrangian equations of motion (|15[) take the following form 

g = - e +n ) + ^ +n (i J+n - f j+n (0) = 0, (22) 
^ = a i+nj (e - e +n ) + a t+n3+n {e +n - f J+n (0) = 0. (23) 

The variables £* +n are auxiliary variables iff <Tj+ n:; + n = 0. We can 
define the second order Lagrangian for the reduced system by apply- 
ing the reduction procedure outlined in section 2. In the case where 
Oi+nj+n ^ the variables £*+ n are not auxiliary and the reduction 
process can not help to construct a Lagrangian for the second order 
system. We will see below how to use the change of variables (|20|) to 
construct an equivalent Lagrangian associated to the reduced origi- 
nal system of second order equations of motion. The two conditions 
upon the transformation (|2(Jj) are 1) <7j-|_ n j-|_ ra (£) = and 2) that the 
transformation be symmetry of the equations of motion. 

Defining the configuration space by the variables £ l we can con- 
struct the explicit relation between the equations of motion associated 
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to the first order system and the reduced space equations (fT2|) . It is 
easy to write this relation using the matrix a. The result is 



(24) 

where cr 4 -? is the inverse matrix of the matrix Oij . This relation suggest 
that the combinations of the Lagrangian equations of motion 



and 



'§ +° ij+n s§^ = e-e +n > ( 26 ) 



are, in fact, equivalent to the system 

f _ = 0, C +n - f +n (£) = 0, (27) 

iff det cr^O. 

But this equivalence is not "variationally admissible" ,2Uj because 
does not define auxiliary variables in the reduction process. This imply 
that the Lagrangian for the second order system ( "Newton equations" ) 
can not be defined by using the first order Lagrangian. 

Notice that here £ l and £ l+n do not have the meaning of coordi- 
nates and momenta, respectively. The nessesary conditions for the 
existence of a first order Lagrangian with auxiliary variables £ 4+n are 
the condition (|17j) in terms of the new variables ^ and a i+n j +n (^) = 0. 

For example consider the isotropic harmonic oscillator in the plane 
with mass and frequency equal to one. The phase space is defined by 
qi, q2,Pi,P2 with the standard symplectic structure. Now consider the 
first order Lagrangian ^!|] 

L = E(qipi - piXii) — H (28) 

where E = \ + Pi) an d H = -2E 2 . The matrix a a b is 

a ( PiQj ~ QiPj 2ES ij + Hlj + PiPj\ (29) 

ab \-(2E5ij +qiqj +piPj) PiQj ~ QiPj )' 

The momenta pi are not more auxiliary variables and the reduction 
process fails. So we can not define the reduced second order La- 
grangian associated to the first order Lagrangian (|28[) by the process 
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of elimination of the variable pi by using its own equations of motion. 
Nevertheless, the dynamical system described by ()28(l is equivalent 
to the original harmonic oscillator. Notice that the phase space in 
this description is highly noncommutative and that a second order 
Lagrangian for the equivalent dynamics, i.e., the harmonic oscillator, 
of course, can be constructed. Another examples along these lines can 
be found in |21j . 

We can use Darboux transformations to construct a Lagrangian for 
a given set of second order equations. The procedure is: 1) Choose the 
variables that we want to eliminate, 2) Construct a Darboux trans- 
formation that renders those variables to auxiliary variables, 3) Elim- 
inate those variables from the first order system and the Lagrangian. 
At the end we will have a Lagrangian for a system of second order 
equations of motion. If this equations are equivalent to the original 
ones we have also an equivalent Lagrangian for the given equations. 
The requirement upon the Darboux transformation is that its projec- 
tion to configuration space (that results in a function of coordinates 
and its derivatives) be a symmetry of the equations of motion, i.e., 
M l = A* M J where M l denotes the original set of equations of motion 
and M? the transformed one and A 1 - a regular matrix. 



In this section we will analyze the case when we have at hand a Hamil- 
tonian function to describe the given system of first order differential 
equations. In this case as we noticed, the matrix (|18j) corresponds to 
the symplectic structure of the Hamiltonian system. So we will ana- 
lyze the role played by the matrix in the definition of the Poisson 
bracket in the associated "phase space" . In particular, we will see that 
the matrix Uij is related with the noncommutativity for the reduced 
coordinates ( "configuration space" ) . In the case where aij / Owe will 
talk about a "noncommutative configuration space" . For that end we 
will split the variables x a into x l ,pj. The matrix a a b can be written 
in the form 



where Bij , dj are antisymmetic matrices that depend on the form 
of the specific Lagrangian that we choose to describe the given first 
order system. The matrix a a b is by construction antisymmetric, reg- 



3 Hamiltonian Systems 




(30) 
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ular det a ^ and closed (a = d£, da = where d is the exterior 
derivative). So it can be used to construct a Poisson bracket 



{F(X) , GW) = _^_ ,31) 

where a ab is the inverse of the matrix a a \,. By construction this Possion 
bracket satisfies the Jacobi Identity (because a = d£). The inverse 
matrix a ab has the form 



-i = ((AT)- 1 CM' 1 A' 1 N~ 1 (A T )~ 1 
° " V -M~ l A- 1 A- 1 BN-\A T )- 1 



(32) 



where M = l+A- l B(A T ) 1 C and iV = 1+(A T ) 1 CA~ 1 B. From here 
we can see the role played by the matrix Cy in the noncommutativity 
of the reduced coordinates x % . Indeed 

{x\x j } = ((/j^CM- 1 ^ 1 )^ (33) 

In this case we say that we have a noncommutative configuration space. 
In particular, when the matrix A = 1 and B = 0, this relation becomes 

{x i ,x j } = C ij . (34) 

For a Hamiltonian system the first order Lagrangian takes the 
general form 

L = £ a (x)x a - H(x), (35) 
and the equations of motion are 

In particular when the Hamiltonian is of the form H = \p 2 + V(x), 
the equations of motion can be written as 



j-i-BaV A l3 i> 3 — v (37) 



SL 

The case C = is 



= (A%*+Cyp j -p j . (38) 



I -A" 1 A- l B{A T ) 1 y J 
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and correspond, in view of the analysis of the previous section, to 
the case where the second order Lagrangian can be constructed be- 
cause the momenta pi are auxiliary variables in the usual sense (usual 
configuration space). In this case we recover the standard reduction 
procedure, even if A, B being arbitrary functions of x 1 . The case C ^ 
is more interesting. In this case, the relation analogous to (|12|) reads 



ah d 
j^j I 

5xi dt 



5L 

dp, 



+ C ik A 



ki 



5L_ 

5x l 



C% k A 



ki 



dV 

dx l 



d r 
+ ~dt 



{{A T ) a + C lk A^B jt )x l 



These equations suggest that the combination 

0, 



— + C ik A kl — 



6p 



5x l 



(41) 



play the role of the "definition of momenta" for this set of equations. 
In fact, the original system of equations is equivalent to the system 



r)V 

C ik A^B^ - C ik A kl — { + s,l'),/ 



r-J 



Bij X^ 



Aypj 



-Pi = 0, 

dV 

dx l 



0, 



(42) 
(43) 



iff det 4/0. By using ()42j) as the definition of momenta we can 
recover the second order dynamics by pluggin this definition in the 
equation (J23J). 

The reduction is, as in the previous case, not "variationaly ad- 
misible" so the reduced Lagrangian can not be constructed by using 
auxiliary variables. 

The second order dynamics is described by the equations of motion 



A* B jk x k 



dxi J 



d_ 

dt 



CikA 



kl 



dV 

dx l 



+ 



dt 



{{A T ) a + C ik A^B fl )x 



(44) 

These "Newtonian equations of motion" must be taken with care be- 
cause the reduction procedure used to obtain them was not standard. 
Indeed, a solution to the original first order system with prescribed 
initial conditions, say by fixing the coordinates and momenta at some 
initial time by x l (to) = Xg,Pj(to) = Pjo, correspond to a solution of 
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the "Newtonian equations of motion" (|44[) with the initial conditions 

x'ito) = 4 (45) 

r 8V i 

C ik A kl B ljX i - CikA^—^ + {A T ) ij V\ (t ) = Pi0 . (46) 

An interesting observation is that these initial conditions are not the 
usual initial position and velocities of Newtonian Mechanics but de- 
pend on the symplectic matrix a and the dynamics through the po- 
tential function V(x). Of course, we can recover the standard case 
when the matrices A = 1, C = 0. 

4 A Non local Lagrangian for the re- 
duced system 

In this section we will construct a non local Lagrangian associated with 
the reduced system (|44|) . Notice that we can not construct a standard 
Lagrangian for this system by the auxiliary variables procedure of our 
previous section. We will restrict ourselves to the case when A,B,C 
are constant matrices. A modification of our procedure in the general 
case is still an open problem. 

We will show here that in this case, we can construct the equations 
of motion in the reduced space and its associated Lagrangian. The 
Lagrangian will be a non local function of the reduced variables and 
its derivatives with respect to time of any order. 

Observing that we can solve in a formal way the equations of mo- 
tion (|38j) for the "auxiliary" variables pj we have 

Pj = (1 " C^Y k (A T )kli l , (47) 

where (1 — C-^y k denotes the inverse of the matrix of (1 — C-^)jk- 
Notice that we have obtained the variables pj that are the analogous of 
the momenta in the standard Hamiltonian formulation in terms of the 
reduced coordinates x 3 and its derivatives with respect to time of infi- 
nite order. This relations can be interpreted in the space of solutions 
as the Green functions associated to the differential equations (|38|) 
given the solution for x 3 (t) with some specified boundary conditions. 
It is possible to extend in some cases the concept of auxiliary variables 
when we have such differential equations and not the algebraic usual 
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ones as @. As far as we know there is no systematic approach for 
this case. We will follow here another route. We will take the for- 
mal solution (|47|) as given and claim that they can be considered as 
"auxiliary variables" . 

With this observations in mind we can use (in a formal way) the 
reduction procedure of section 2. In particular the equations of motion 
in the reduced space are 



A d 
Atk Jt 



and by multiplication of the matrix (1 — C^)jk this system is equiv- 
alent to the system (|44|). 

To construct the reduced Lagrangian function we observe that the 
first order Lagrangian associated to the system (|37I38|) is 

L = ^a ab x a x b -H(x), (49) 

where x a = (x l ,pj) and the matrix a ab is given by (|3U|) , In terms of 
x l ,pj this Lagrangian is 

L = ^B ij x i x j + AijPji* + ^CijPipi - H. (50) 

Pluggin pj from (|47|1 into this Lagrangian function we obtain the re- 
duced non local Lagrangian 



L NL = - 2 lhjy-r> + \Aijx\l - C^y k A T kl x l - V(x). (51) 



This reduced Lagrangian function is also of infinite order in time 
derivatives of the reduced coordinates x % and to obtain from it the 
equations of motion Q48|) we need to apply the extended Euler-Lagrange 
operator 

5 — = V( 1)^ W dL (52) 

where (x 1 )^ denotes the n— derivative with respect to time of the 
coordinate x % . 

The equations of motion of the Lagrangian (|51|) using the opera- 
tor (|52|) are precisely l|48|). The second order dynamics (|44|) is equiv- 
alent to the equations of motion Q48JI associated with the non-local 
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Lagrangian (|51j). i.e., 



fl - C— 



d_ 
dt 



CikA' 



hi 



dV 
dx l 



d_ 

+ Jt 



((A^a + Qk^fl,,)* 1 



(53) 



As the matrix operator (1 — C4r)ij have a formal inverse (to obtain 
the momenta 1)47(1 ) the system f)53f) is equivalent with the system (|44j) . 



5 Noncommutative Classical Dynam- 
ics in configuration space 

In this section we will develop, as a simple application of the ideas 
above, the noncommutative classical mechanics constructed, as we 
stated in the introduction from a classical Hamiltonian function and a 
constant symplectic structure. We will outline some properties of this 
noncommutative dynamics and the reduced non-local Lagrangian. 

Lets start by the definition of the noncommutative classical me- 
chanics as the dynamics described by the first order Lagrangian 



yp y + 0p y p x - axy) - H(x, y,p x ,P y ) 



(54) 



L = a (xp x 

where a = 1 } ct0 . In the case when H = \p 2 + V(x,y) the associated 
Lagrangian equations of motion are 



5L 

- — = a(x + 

OPx 



SL 

bp y 



) ~Px 

a{y - Qpx) ~Py 



5L . . 

— = a{-p x + ay) 
Ox 

SL 

— = a(-py - ax) 
by 



dV 

dx 
dV 

dy 



0, 
0. 
0, 
0. 



(55) 
(56) 

(57) 
(58) 



So the matrix a is 



a 



a. 



( a -1 \ 

-a -1 

1 6* 

\ 1-61 / 
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from our previous notation we can identify 

A = a G i)' B = a (-° o)' c = a {-e o 

The relation ()12j) is in this case 

d /<5L\ /<JL\ /<5£\ „ ay+ ldV e d ( dV ) 

dt \ 5p x J dt \Sy J \Sx J a dx dt dy 

d{5L\ n d fSL\ {5L\ 1 dV n d,dV, 



dt \ 5py J dt \5x J \Sy J a dy dt dx 

From here we observe that there is a combination of equations of 
motion in "phase space" that can play the role of the "definition of 
momenta" . The combination is 

SL JL , . 

(59) 

-£-#-0. (60) 

OPy OX 

Using these relations the system 

1 dV 

fc-^ + -^ = 0, (61) 
1 dV 

p y + ax + -— = 0, (62) 
a ay 

dV 

p x = x-B—, (63) 
dV 

Px = ij + e ~dy~ 1 (64) 

is equivalent to the original system. The key point to notice here is 
that the definition of the momenta is not purely kinematics but con- 
tains dynamical information, i.e., are different for each system defined 
by some potential function V(x, y). This implies that the second order 
dynamics given by 

ldV n d,dV, , , 

+ -»*<«»> = °- (65) 

ldV n d,dV, . . 

! ' +ra + ^ + V&) = ' < 66 » 

is not directly equivalent with the original first order system. We need 
to specify appropiarte boudary conditions to associate to a given so- 
lution of the original first order system the corresponding solution for 
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the second order system (|65l New-eq-y). For example the Lagrangian 
1)541) implies the following fixing at the boundaries 

9 9 

(x + ^Py)(U) = Xi, (x + ~Py)(tf) = Xf, (67) 

9 

{y - ^Px)(U) = yu (y - -p x )(t f ) = y f . (68) 

For each trajectory with the specified boundary conditions ()67)1 and 
(|68j) corresponds a solution of the system 1)65)1 and ()66)) with the 
boundaries 

9 dV \ / 9 dV \ 

+ 2 (V + ^)) = [ x + 2 ( ^ + ^) = */> ( 69 ) 

(v-2 (i_ ^ ) ) ( * i) = W ' ( y -2 (i -W) ( * /)=y/ - (70) 
This is the exact sense in which the original and the reduced second 
order dynamics are equivalent. It is worth noting that these boundary 
conditions depends on the parameter 9 and on the specific potential 
that appears in the Newtonian equations (|65j) and ()66j) . This infor- 
mation is crucial for the quatization of the classical system using the 
Path Integral. 

To construct the non-Local Lagrangian associated to the equations 
of motion 1)65)1 and ()66j) we obtain the momenta from equations 1)55)1 
and ()56)) to get 

d 2 _ d d 2 _ d 

p x = a(l+(a9) 2 -^)~ 1 — (x+a9y), p y = a(l+(a9) 2 — (y-a9x) 

(71) 

According to our analysis the equations of motion in the noncom- 
mutative configuration space are 

hit ~ i* + (1 + {ae) 'w' rl ^ {x + a6i) = (72) 

and the Lagrangian is 

L = ^(x*x + y*y) + ^(x*y^y*x) + ^(xy-yx)-V(x,y), (74) 
where 

iSQ^l + JaS) 2 ^)" 1 . (75) 

The equations ()72j) and l)73|) are equivalent to the second order 
equations (jfjjj) and ()66)) . 
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6 Examples 



Consider the following Lagrangian 16 , 

L = (y + z)x + wz + — (w 2 — 2zy — z 2 ). 

The symplectic structure that can be read out from this Lagrangian 
is 

{x,y} = l, {y,w} = -l, {w,z} = -l, 
(all other Poisson brackets are zero). The equations of motion are 

g = -V-i = 0, (76) 
£-*-«-<>, (77) 

A T 

= z + w = 0, (78) 

ow 

S L 

— = x — w — y — z = 0. (79) 

oz 

At first sight the variables y, z as a set and w are auxiliary. In fact, 
using (|77 j) .(|79 |) and (|78|) we can define a reduced Lagrangian for the 
variable x that gives the correct reduced Lagrangian and four order 
equation of motion for x. We can also consider as auxiliary vari- 
ables only y, z leaving the reduced space for x, w whose Lagrangian 
can also be constructed. In this case we observe that the reduced 
coordinates define a commutative configuration space (x, w have zero 
Poisson bracket). 

Now suppose that we want to construct the reduced dynamics 
in the noncommutative space x,y. To arrive to this space we can 
use the equations (f77|) and (|T5|) to eliminate the variables z, w but 
this reduction is not variational admissible so we can not construct 
a Lagrangian for this reduced system. In fact, it can be shown that 
a quadratic Lagrangian for this reduced second order dynamics does 
not exist. 

The construction of the non-local Lagrangian can be worked out 
by the procedure sketched in this note and consist in the use of the 
equations of motion (|79|) and (|78|) to obtain the variables z, w from its 
own equations of motion. The first observation is that these variables 
are not auxiliary in the usual sense, but nevertheless we can obtain 
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them as a formal integration of the equations (|79[) and (|78[) to get 

*=p-^)-'(*-»). ---d-^)- 1 ^*-*)- (») 

The reduced equations of motion are 

i - (i - - v) = °> y + (i - ^y l j t ^ ~ y) = °- ( 81 ) 

These equations are equivalent to the second order equations 

y + y = 0, x + y = 0. (82) 

The reduction process, namely the construction of a Lagrangian asso- 
ciated to this equations of motion can be implemented to get 

L R = ^(±-y)(l-^r 1 (x-y) + y±. (83) 

This Lagrangian reproduces the equations of motion (|81|). The re- 
lation between the Lagrangian equations of motion and the reduced 
second order dynamics is given by 

d 5L 5L , 
dtJy- + T X = X + y ' (84) 
d 5L d 5L SL + 6L _ + 
dt 5x dt 5w bz by 

From here we can see that the combination |^ + = and ^ = 
define the "momenta" 

w — y = 0, x — z = 0. (86) 

The system 

w — y = 0, x — z = 0, z + w = 0, x — w — y — z = 0, (87) 

is equivalent to the original set of first order equations of motion. In 
this case the general solution of the original first order system and 
the general solution for the reduced system Q82JI coincides because the 
"momenta" z, w as given by (|8(ij) are precisely the velocities as in the 
standard Newtonian case. 

The case of the noncommutative classical mechanics analyzed in 
section 5 can be illustrated by taking as an example the simple poten- 
tial V = c{x+y) with c any constant. Here the general solutions of the 
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original first order system and the reduced one do not coincide, but 
taking into account the properly choosed boundary conditions ()69|) 
and Q7UJ) the respective solution coincides and are 

x{t) = — ^-4tfXi + 4tiXf + AtX -29Y - 
cT(2t 2 + 0(t f + 0)- 2t(U +t f + 9) + ti(2t f + l 



cT[2t 2 + 6(6-t f )- 2t(U + t f -6) + ti(2t f 



y(t) = — ( y -4t f y l + U i y f + 4tY + 29X 



where X = x% — Xf, Y = yi — yf, T = ti — tf. The free particle can be 
obtained by taking c = in the expressions above. It is interesting to 
observe that due to the boundary conditions the solution depends on 
the parameter 9. 

As a final example consider the construction of a second order 
Lagrangian for the noncommutative harmonic osillator described by 
the first order Lagrangian (|54|l. with the potential given by V = 
^ui 2 (x 2 +y 2 ) , using a Darboux transformation to render some variables 
to auxiliary variables. The Darboux transformation is 

x' = x + ^9p y , y' = y-^9p x , p' x =p x , p' y =Py (88) 

With this transformation we construct a new first order Lagrangian L' . 
For this Lagrangian the momenta p' x ,p' y are auxiliary variables. Then 
we proceed to eliminate them to find the second order Lagrangian 



L' = n 2 



2 



I -12 , -/2\ 1 2/ 12 , I2\ i ^ / 

(x +y {x +y ) + —{xy -yx) 



where = (1 + (^) 2 ) -1 - The crucial point is that the equations of 
motion associated to this second order Lagrangian are equivalet to 
the original second order system. The equations of motion associated 
with this Lagrangian are 

n 2 [x' + u?x' + uj 2 9y'] = 0, 

n 2 [y' + u?y + lo 2 9x'] = 0. 

These equations of motion are proportional to the original system of 
second order equations. So the Darboux transformation (|88|) is a sym- 
metry of the equations of motion. So, we have constructed an equiva- 
lent Lagrangian for the given system of second order equations. The 
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symmetry transformation is given by the projection of the Darboux 
transformation to the configuration space 

x = x'- \etf{y> + ^-x'), y = y>- \etf(x' + ^y'). 

7 Conclusions 

In summary, we have analyzed the role of auxiliary variables in first or- 
der dynamics and its relation with the inverse problem of the calculus 
of variations. The relation between the first order equations and the 
reduced second order equations was worked out by showing that there 
exist a relation between a functional combination of the first order 
equations of motion and the reduced system. This relation can always 
be computed when the variables of the reduced space are fixed and 
works a the level of the equations of motion. When the eliminated 
variables are auxiliary, the Lagrangian can be constructed directly. 
When the variables are not auxiliary we need to find a symmetry of 
the equations of motion that allows the construction of a Lagrangian 
with auxiliary variables. In this way we can construct an equivalent 
second order Lagrangian for a given set of second order equations of 
motion. 

When the symplectic structure is constant we can construct, in 
a formal way, a non-local Lagrangian. The equations of motion of 
this non-local Lagrangian are equivalent to the set of second order 
equations (|44[). These equations are in turn equivalent to the original 
set of first order equations of motion. 

Finally we have worked out as examples the case of "non-commutaive 
mechanics" , a system with no second order Lagrangian and for the case 
of the non-commutaive harmonic oscillator we have constructed using 
our procedure an equivalent Lagrangian for the original second order 
dynamics. It will be interesting to analyze the problem of the quanti- 
zation of this systems and its relation with the non-commutaive 
quantum mechanics. 
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